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Summary 

The African clawed frog, Xenopus laevis, is an ectothermic 
vertebrate that can survive at low environmental 
temperatures. To gain insight into the molecular events 
induced by low body temperature, liver proteins were 
evaluated at the standard laboratory rearing temperature 
(22 °C, control) and a low environmental temperature (5°C, 
cold exposure). Using nano-flow liquid chromatography 
coupled with tandem mass spectrometry, we identified 58 
proteins that differed in abundance. A subsequent Gene 
Ontology analysis revealed that the tyrosine and 
phenylalanine catabolic processes were modulated by cold 
exposure, which resulted in decreases in hepatic tyrosine and 
phenylalanine, respectively. Similarly, levels of pyruvate 
kinase and enolase, which are involved in glycolysis and 
glycogen synthesis, were also decreased, whereas levels of 
glycogen phosphorylase, which participates in glycogenolysis, 
were increased. Therefore, we measured metabolites in the 
respective pathways and found that levels of hepatic glycogen 



and glucose were decreased. Although the liver was under 
oxidative stress because of iron accumulation caused by 
hepatic erythrocyte destruction, the hepatic NADPH/NADP 
ratio was not changed. Thus, glycogen is probably utilized 
mainly for NADPH supply rather than for energy or glucose 
production. In conclusion, X. laevis responds to low body 
temperature by modulating its hepatic proteome, which 
results in altered carbohydrate metabolism. 

© 2013. Published by The Company of Biologists Ltd. This is an 
Open Access article distributed under the terms of the Creative 
Commons Attribution License (http://creativecommons.org/ 
licenses/by /3.0), which permits unrestricted use, distribution 
and reproduction in any medium provided that the original 
work is properly attributed. 
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Introduction 

In order to survive, organisms respond to environmental changes 
by altering their expression patterns of genes and proteins. One of 
the most important environmental factors is temperature, which 
governs the distribution, behaviour, and physiological response 
of organisms. When the environmental temperature falls, 
endothermic vertebrates, i.e. mammals and birds, maintain their 
core body temperature within a narrow range by a 
thermoregulatory system. Their bodies undergo numerous 
physiological changes in response to cold stress to maintain 
their temperature (Sonna et al., 2002; Silva, 2006). A cold 
environment leads to hypothermia - abnormally low body 
temperatures. Apart from species that hibernate, most 
endotherms are unable to survive at low body temperatures. In 
contrast to endotherms, the core body temperature of ectothermic 
vertebrates, including most fish, amphibians, and reptiles, as well 
as most invertebrates, becomes low (Salt, 1949). The low body 
temperature conditions modulate their physiological functioning 



directly or indirectly due to low environmental temperature. How 
and why mammalian hibernators and ectothermic vertebrates 
survive in low temperatures have long been a scientific challenge 
and an important line of inquiry for biologists. The cellular and 
physiological responses to low body temperature have been 
investigated. The key response is the downregulation of the 
cellular metabolic rate to new hypometabolic steady states in a 
way that balances the ATP demand and ATP supply pathways 
(Hochachka, 1986; Boutilier, 2001). 

Amphibians have been used as physiological models to study 
the responses to environmental stresses (Burggren and 
Warburton, 2007; Hopkins, 2007). Cold exposure depresses 
their metabolic rates and has effects on mitochondrial 
bioenergetics (Boutilier et al., 1997; Trzcionka et al., 2008). 
The African clawed frog, Xenopus laevis (Daudin), has served as 
a crucial vertebrate model for biological research. Since the wild 
habitat distribution of X. laevis covers a wide geographical range 
(Tinsley et al., 1996; Tinsley and McCoid, 1996), this species is 
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capable of tolerating a range of environmental conditions. 
Therefore, we considered that X. laevis exposed to low 
temperatures would be a good hypothermic model to 
investigate responses to environmental temperature. 

To understand the physiological responses to low-temperature 
stress, global analyses of the genes and proteins responding to an 
environmental stimulus should enable a deciphering of the 
underlying physiological pathways (Cossins et al., 2006; Gracey, 
2007). Although a number of previous studies have compared the 
expression levels of mRNAs and proteins in cells and tissues, their 
correlation has been controversial (Anderson and Seilhamer, 1997; 
Ideker et al, 2001; Chen et al, 2002; Griffin et al, 2002; Mehra et 
al, 2003; de Sousa Abreu et al, 2009; Maier et al, 2009; 
Schwanhausser et al, 2011; Ghazalpour et al, 2011). Considering 
that physiological events are determined by protein-driven 
processes, the proteome should provide key information to 
understand molecular responses. Proteomic studies have 
investigated the following associated physiological responses in 
X. laevis: the biosynthetic and secretory processes in 
neuroendocrine cells (van Herp et al, 2008), white/black- 
background adaptation (Devreese et al, 2010), response to 
chemical exposure (Gillardin et al, 2009; Serrano et al, 2010), 
photoreceptor outer segment assembly (Wang et al, 2009), and 
limb regeneration (King et al, 2009). DNA microarrays forX. laevis 
studies are currently commercially available from some companies 
(GeneChip Xenopus laevis Genome Array from Affymetrix, Inc. 
and Xenopus Gene Expression Microarray from Agilent 
Technologies, Inc.). However, the proteomic approach mX. laevis 
has been limited due to lack of availability of complete genomic 
information. The draft genome sequence assembly of the Western 
clawed frog, Xenopus (Silurana) tropicalis, was recently reported 
(Hellsten et al, 2010). Additionally, concerted efforts are currently 
underway in the United States, Japan, and elsewhere to carry out X. 
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laevis genome sequencing [Xenopus Community White Paper 
2011; Xenbase: Xenopus laevis and Xenopus tropicalis biology and 
genomics resource (http://www.xenbase.org/common)]. In the near 
future, improvement in the accuracy of the Xenopus genome 
sequence will be achieved to facilitate proteomic approaches. 

We previously investigated the haematopoietic response to low 
temperature in X. laevis because haematopoiesis is one of the 
most important physiological functions. After 24 hours of cold 
exposure (5°C), X. laevis displays anaemia associated with 
hepatic erythrocyte destruction and hepatic iron accumulation as 
a result of heme degradation (Maekawa et al, 2012). The 
anaemia is prolonged during cold exposure concomitantly with 
hepatic confinement of newly produced erythrocytes (Maekawa 
et al, 2012). Generally, the liver plays a central role in metabolic 
homeostasis and is a major site for the synthesis, metabolism, 
storage, and redistribution of carbohydrates, proteins, and lipids 
(Bechmann et al, 2012). The liver also plays an important role in 
energy metabolism and the large change in metabolic rate caused 
by cold exposure. In X. laevis, in contrast to mammals and 
terrestrial frogs, the liver also plays a central role in the 
production and destruction of erythrocytes (Chegini et al, 1979; 
Nogawa-Kosaka et al, 2010; Nogawa-Kosaka et al, 2011; 
Maekawa et al, 2012; Okui et al, 2013). Therefore, after 
24 hours of cold exposure, various physiological responses may 
occur in the liver. It is important to investigate these initial 
responses to cold-exposure-induced low body temperature. 

In the present study, we used a proteomics approach to profile 
the liver proteome in X. laevis after exposure to low temperature, 
because proteomics studies on the liver at low body temperature, 
such as gilthead sea breams exposed to the cold (Ibarz et al, 
2010), mammalian hibernators during entrance into hibernation 
(Epperson et al, 2004; Shao et al, 2010; Epperson et al, 2010; 
Rose et al, 2011), rat induced hypothermia (Oda et al, 2012), 
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Fig. 1. Transitions in environmental temperature and 
hepatic proteins. (A) Relationship between cold-exposure time 
and water temperature. A plastic tank containing 1 1 of 22 °C 
water was transferred to an incubator at 5°C and allowed to cool. 
During that time, the temperature of the water was measured. (B) 
Image of SDS-PAGE gels used to separate X. laevis liver protein 
samples; 1 0 u.g of each sample from either the control or cold- 
exposure group (three frogs each) was separated on 12% SDS 
polyacrylamide gels. Protein bands were visualized by staining 
with Coomassie brilliant blue R-250. Reducing (right panel) and 
non-reducing (left panel) conditions are shown. 
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and freeze-tolerant wood frogs during winter (Kiss et al., 2011) 
were recently reported. These studies enable us to carry out cross- 
species comparisons of liver proteome changes. We applied a 
label-free quantification method using nano-flow liquid 
chromatography coupled with tandem mass spectrometry 
(nanoLC-MS/MS) to assess X. laevis liver proteins that differ 
in abundance between standard laboratory rearing temperature 
(22 °C, control condition) and low environmental temperature 
(5°C, cold exposure). The purpose of this study was to gain an 
insight into the initial physiological response to cold-exposure- 
induced low body temperature. 

Materials and Methods 

Animals 

Wild-type X. laevis (mass 30^0 g) frogs were purchased from Kazuo Ouchi 
(Misato, Saitama, Japan) and housed in plastic tanks at the standard laboratory 
rearing temperature (22 °C) with constantly running water. This condition was 
defined as the control condition. For low-temperature exposure, plastic tanks 
containing X. laevis (one frog per tank containing 1 1 of water at 22 °C) were 
transferred to an incubator (Bio Multi incubator; NK Systems, Osaka, Japan) set at 
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5°C and allowed to cool. All experiments were conducted according to the 
Regulations for Animal Experimentation at Waseda University. 

Haematological analysis 

The haematological value of peripheral blood, including blood cell counts, 
haemoglobin, and haematocrit values were obtained as previously reported 
(Aizawa et al., 2005; Nogawa-Kosaka et al., 2010; Nogawa-Kosaka et al., 2011; 
Maekawa et al., 2012). 

Liver tissue collection 

At 24 hours after cold exposure, X. laevis were quickly killed by beheading. The 
livers from each of the control and cold-exposure groups (#=3 each) were cut into 
smaller pieces and flushed with Tris-buffered saline (20 mmol l _l Tris-HCl 
pH 7.5, 100 mmol 1 _1 NaCl) containing 1 mmol 1 _1 ethylenediaminetetraacetic 
acid (TBSE) to remove the excess blood, and then directly frozen in liquid nitrogen 
and stored at — 80°C for protein analyses. 

Extraction of hepatic proteins 

A piece of the liver was homogenized in TBSE (0.5 ml per 100 mg of tissue) using 
a bead beater-type homogenizer (Beads Crusher uT-12; TAITEC Co., Ltd., 
Saitama, Japan). The liver homogenates were centrifuged at l,500xg for 5 minutes 
at 4°C to remove cell debris and the supernatants were further centrifuged at 
15,000xg for 20 minutes at 4°C to remove insoluble proteins. The supernatants 
containing soluble protein were collected and stored at — 80 °C until use. Protein 
concentration was determined using the Bradford assay reagent. 

Protein digestion 

Three protein extracts from each individual were mixed in equal amounts 
(Fig. 2). The mixed extract containing 50 jig of protein was dissolved in 
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Fig. 2. Schematic view of sample and data processing. Each test sample, i.e. 
a protein extract mixture derived from three frogs prepared by bead crusher, 
was subjected to tryptic digestion. Subsequent nanoLC-MS/MS analysis was 
repeated three times. Automatic peak extraction and MASCOT MS/MS ion 
search were performed using Xome software. The triplicate data were 
processed individually in all combinations to calculate the protein ratios by 
non-label quantification using Mass Navigator software then protein ratios were 
averaged after normalization. The details are explained in the Materials and 
Methods and in Table 1. 



Fig. 3. Outlines of MASCOT MS/MS ions search. (A) Summary of 
MASCOT peptide identification in each search. "Peptide matches above 
identity threshold (P<0.05); *FDR = decoy hits/NCBInr hits. (B) Venn 
diagrams of identified proteins in individual searches. Proteins identified in at 
least two of the triplicate runs are highlighted in grey. (C) Comparison of valid 
identified proteins differentially expressed between the control and cold- 
exposure conditions. Proteins identified at least twice were compared (126 
proteins from the control group and 100 proteins from the cold-exposure 
group); 8 1 proteins overlapped between both groups. 
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0.5 mol 1 1 Tns-HCl (pH 8.5) containing 8 mol 1 1 urea, 2.5 mmol 1 1 
ethylenediaminetetraacetic acid, and 10 mmol l _l dithiothreitol, and incubated for 
1.5 hours at 37°C. Iodoacetamide was then added at a concentration of 50 mmol 1 _1 to 
alkylate the reduced thiol groups. After incubation for 30 minutes at room temperature 
in the dark, the mixture was diluted with 50 mmol 1 _1 ammonium bicarbonate buffer at 
a final concentration of 1 mol 1 1 urea. For digestion, sequencing grade modified 
trypsin (Promega KK., Tokyo, Japan) was added to the protein solution at 1:50 
(trypsin: protein) and the mixture was incubated for 19 hours at 37°C. Formic acid was 
added at a concentration of 0.1% to stop the reaction. The final solution was cleaned up 
with MonoTip C18 tips (GL Sciences Inc., Tokyo, Japan) and the eluent was 
evaporated in a vacuum centrifuge. The powdered peptides were redissolved in 50 Jll 
of 2% acetonitrile in water containing 0.02% formic acid for nanoLC-MS/MS analysis. 

nanoLC-MS/MS analysis 

We used Nano Frontier eLD system (Hitachi High-Technologies Corporation, 
Tokyo, Japan). The analytical column was a packed nano-capillary column 
(NTCC-360/75-3; Nikkyo Technos Co. Ltd., Tokyo, Japan). We also used a 
monolithic trap column (Monolith Trap C 1 8-50- 15, Hitachi High-technologies 



Corporation, Tokyo, Japan). The LC conditions were as follows. The flow rate 
of the nanoflow pump was set at 200 nl min — , Solvent A was 2% aqueous 
acetonitrile containing 0.1% formic acid, and solvent B was 98% aqueous 
acetonitrile containing 0.1% formic acid. The composition of solvent B was 
linearly increased from 2% at 0 minutes to 35% at 150 minutes, maintained at 
100% until 165 minutes, and then returned to the initial condition of 2%. The 
MS/MS conditions were as follows: ESI voltage, +1.6 kV; curtain (nitrogen) 
gas flow rate, 0.6 1 min - ; precursor mass scan range, m/z 100-2000; scan 
time, 20 msec; fragment mass scan range, m/z 50-2000. The 2 (il of peptide 
solution (equivalent to 2 fig of protein) were analysed in triplicate runs 
(Fig. 2). 

Data processing 

Xome (Mitsui Knowledge Industry Co., Ltd., Tokyo, Japan) software was used for 
generating the peak list and identification of proteins (Honmyo, 2007). For the 
identification of proteins, we performed peptide mass fingerprint against the 
NCBInr database [NCBInr 20130303 (23463169 sequences; 8064228071 
residues); National Center for Biotechnology Information, http://www.ncbi.nlm. 



Table 1. Upregulated proteins in cold-exposed X. laevis liver (group 1). 
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a Protein scores and peptide contents in MASCOT MS/MS ion search. Maximum values for triplicate data are shown. 

b Fold changes (cold exposure/control) of protein abundance. Identified peptides were aligned by m/z and retention time. Peptide peaks of the MS chromatogram 
were approximated by the Gaussian waveform using the Simplex method and peak areas were calculated. Peptide ratios (cold exposure/control) were calculated 
from the peak areas. Outliers were removed by Thompson's rejection test and then protein ratios were estimated as the average of the peptide ratio. The triplicate 
data were processed individually in all combinations to calculate protein ratios, and then the protein ratios were averaged after normalization. In many proteomics 
studies, the expression of proteins does not change under different conditions (Tabata et al., 2007), and the mode value (i.e. the most probable value) of protein 
ratios reflects variations in protein loading. Therefore, to take the variations into consideration, the ratio (cold exposure/control) for individual proteins between 
samples was normalized against the mode value of the ratios (the mode value was set as 1.0 for each combination). The scheme of this procedure is described in 
Fig. 2. 

By reference to a previous publication (Kolkman et al., 2006), the ratios were considered as significantly changed when they satisfied the following assumption: 
proteins with a cold-exposure/control ratio higher than 1.25 or lower than 0.8 had to have a coefficient of variation, i.e. average ratio divided by standard 
deviation, below 30%. Significant values are in bold. The information for the peptide used for calculating the protein fold change (e.g. retention time and peak 
area) is provided as supplementary material Table S5. 

c RefSeq IDs of human homologues were obtained from NCBI HomoloGene database or by alignment with the NCBI Reference Sequence (RefSeq) database 
using the BLASTp programme (^). 
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nih.gov] using the Mascot search engine version 2.1 (Matrix Science, London, 
UK). The Mascot MS/MS ion search conditions were as follows: taxonomy filter, 
X, laevis (African clawed frog) (17403 sequences); enzyme, trypsin; maximum 
missed cleavages, I; fixed modifications, carbamidomethyl (C); variable 
modifications, oxidation (M); peptide mass tolerance, 0.3 Da; MS/MS ion mass 
tolerance, 0.3 Da; charge states, +1, +2, and +3; mass values, monoisotopic; 
instrument type, ESI-TRAP. The threshold score to achieve P<0.05 was set by the 
Mascot algorithm. Because the nanoLC-MS/MS was performed in triplicate, the 
triplicate data were processed separately (Fig. 2). As a result, the false discovery 
rate (FDRs) was less than 0.05 in all Mascot searches, which indicates that the 
significance threshold (P<0.05) was applicable (Fig. 3A). The complete results 
are listed in supplementary material Table SI. Only an identification observed in at 
least two of the three replicates was taken to be valid (Fig. 3B). The Xome 
programme identified 145 X. laevis proteins, 113 of which had a human 
homologue in the NCBI HomoloGene database (http://www.ncbi.nlm.nih.gov/ 
homologene). The remaining 32 proteins were manually annotated by alignment 
with the NCBI Reference Sequence (RefSeq) database using the BLASTp 
programme (NCBI Basic Local Alignment Search Tool, http://blast.ncbi.nlm.nih. 
gov/Blast.cgi). 

To screen proteins differentially expressed between the control and cold- 
exposure groups, fold changes in protein abundance were calculated using the 
'non-label quantitation' function of Mass Navigator vl.2 (Mitsui Knowledge 
Industry Co., Ltd., Tokyo, Japan). The procedure for the calculation of protein fold 
changes is described in Table I. The complete results are listed in supplementary 
material Table S2. 



Glycerol Colorimetric Assay Kit (Cayman Chemical Company, MI, USA). Liver 
extracts for the glycerol assay were prepared as previously described (Driedzic et 
al., 2006). A piece of frozen liver was homogenized in nine volumes of 10% 
perchloric acid, and the homogenate was centrifuged at 15,000xg. The supernatant 
was neutralized using potassium hydrate and assayed as described above. The 
absorbance at 540 run was measured. NADP and NADPH in the liver were 
measured using a fluorimetric assay (SensoLyte NADP/NADPH Assay Kit; 
AnaSpec, Inc., CA, USA) according to the manufacturer's instructions. The 
resulting red fluorescence was monitored at excitation and emission wavelengths 
of 560 and 590 nm, respectively. 

Results 

Anaemia is induced by cold exposure in X. laevis 
Over 6 hours of exposure to 5 °C, the water temperature decreased 
to almost 5°C (Fig. 1A). At 24 hours, the erythrocyte count was 
approximately 70% of that in the control condition, and 
haemoglobin and haematocrit values were also decreased as 
reported previously (Maekawa et al., 2012) (data not shown). No 
significant changes in body or hepatic weight were observed 
(data not shown). 



Gene ontology and pathway analysis 

To identify significantly represented biological themes and functional groups in 
the protein list, gene ontology (GO) and pathway analysis were performed using 
the Database for Annotation, Visualization and Integrated Discovery (DAVID) 
v6.7 program (http://david.abcc.ncifcrf.gov) (Huang et al., 2009a; Huang et al., 
2009b). The GO analysis was used to identify enriched biological themes using 
GO terms defined and provided as official terms by the Gene Ontology 
Consortium (http://www.geneontology.org) (Dennis et al., 2003). The 'biological 
process' subontology of GO (GO:BP) refers to a biological objective to which the 
protein contributes and is widely used to evaluate sets of relationships between 
proteins. The pathway analysis was used to identify candidate proteins involved in 
pathways from the Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway 
database (http://www.genome.jp/kegg/pathway.html). The analysis conditions 
were as follows: when the list of identified proteins in the control and cold- 
exposure conditions was used as the input data, the DAVID default population 
background (corresponding genome-wide genes with at least one annotation in the 
analysing categories) was used; when the differentially expressed protein list was 
used, the total protein list from this study (145 proteins) was inputted and used as a 
customized population background. We used EASE scores, which modified 
Fisher's exact test P values to evaluate the significance of enrichment (Hosack et 
al., 2003), and Benjamini and Hochberg FDR procedures for multiple testing 
corrections (Benjamini and Hochberg, 1995). Only statistically enriched GO terms 
and pathways (Benjamini and Hochberg FDR-corrected P<0.01) with at least two 
proteins were selected. To determine the hierarchical structure of the selected GO 
terms, we used QuickGO (http://www.ebi.ac.uk/QuickGO). 

Quantification of free amino acids 

Free amino acids were quantified by EZ:faast GC/FID Free (Physiological) Amino 
Acid Kit (Phenomenex, Inc., CA, USA) by using gas chromatography (GC) 
equipped with a flame ionization detector (FID) (GC-2014; Shimadzu Corporation, 
Kyoto, Japan). The protein extracts from X. laevos liver (equivalent to 0.5 mg 
protein) were subjected to pretreatment for GC/FID analysis according to the 
manufacturer's protocol. Chromatographic conditions were as follows: initial oven 
temperature of 80 °C with 20 °C min~ ramps to the final temperature of 320 °C. 
Inlet temperature of 280 °C under constant helium flow of 50 cm sec -1 and the 
detector temperature of 320 °C is used for the analysis of derivatized amino acids. 

Quantification of glycogen in the liver 

Liver glycogen was isolated by precipitation from ice-cold 5% trichloroacetic acid 
extracts using ice-cold 70% ethanol as previously described (Graff and Allen, 
1963). The precipitate was dissolved in water, hydrolyzed to glucose for 2 hours at 
100 °C in 2 N sulphuric acid, and neutralized using 2 N sodium hydrate (Sahyun, 
1931). The glucose was measured using the glucose oxidase method (Glucose CII 
Test-Wako; Wako Pure Chemical Industries, Ltd., Osaka, Japan). The absorbance 
at 505 nm was measured using a microplate reader (POWERSCAN HT; DS 
Pharma Biomedical Co., Ltd., Osaka, Japan), and a factor of 0.927 was used to 
convert the values obtained from glucose into values for anhydrous glycogen. 

Quantification of glucose, glycerol, NADP and NADPH 

Glucose in the plasma and liver extracts was measured directly using the Glucose 
CII Test from Wako. Glycerol in the plasma was measured directly by using the 
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Fig. 4. GO and pathway analyses ofX. laevis liver proteomes. (A) Schematic 
view of the data processing procedure. The Gl accessions of X. laevis proteins 
were converted to human RefSeq protein IDs by using NCBI HomoloGene and 
BLASTp, and then GO and pathway analyses were performed by using the 
DAVID program. FDR-corrected P values were defined by modified Fisher's 
exact test with the Benjamini and Hochberg FDR correction. The significantly 
identified GO:BP terms appearing deepest in the hierarchy and the significantly 
identified KEGG pathways are shown (FDR- corrected / ? <0.01). The details are 
explained in the materials and methods. (B,C) Comparison of enriched GO:BP 
terms (B) and KEGG pathways (C) in the list of proteins identified between the 
control and cold-exposure conditions. 
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X. laevis liver proteome 

To survey hepatic protein profiles of X. laevis in the control and 
cold-exposure groups, we conducted a proteomics analysis using 
label-free nanoLC-MS/MS. Although the SDS-PAGE patterns of 
the Coomassie brilliant blue-stained liver extracts were similar 
between control and cold-exposure groups (Fig. IB), the LC- 
MS/MS analysis detected differences in the protein contents. We 
identified 145 proteins (126 proteins from the control group and 
100 proteins from the cold-exposure group), and 81 proteins 
overlapped between the groups (Fig. 2C). Next, we set out to find 
biologically relevant themes within the list of identified proteins. 

The GO analysis identified five and seven GO:BP terms that 
were significantly enriched in the control and cold-exposed X. laevis 
liver proteomes, respectively. Four of these terms were common 
between control and cold exposure. Only one term, 'oxygen and 
reactive oxygen species metabolic process (GO:0006800)' (4 
proteins, P=2.87xl0~ 3 ), was control-specific, and three terms, 
'cellular amide metabolic process (GO:0043603)' (4 proteins, 



P=l. 17x10 3 ), 'coenzyme metabolic process (GO:0006732)' (6 
proteins, P=5.89xl0~ 3 ), and 'cellular amino acidmetabolic process 
(GO:0006520)' (9 proteins, P=6.73xl0~ 4 ), were cold-exposure- 
specific. In the pathway analysis, two pathways, 'glycolysis/ 
gluconeogenesis (xlaOOOlO)' (control: 15 proteins, P=2.62xl0~ 9 ; 
cold exposure: 12 proteins, P=3.63x 10~ 7 ) and 'pentose phosphate 
pathway (xla00030)' (control: 9 proteins, P=8.29xl0~ 7 ; cold 
exposure: 9 proteins, P=3.48xl0~ 7 ), were commonly associated 
with the control and cold-exposed liver proteomes mX. laevis. Three 
pathways, 'tyrosine metabolism (xla00350)' (6 proteins, 
P=1.27xl0~ 3 ), 'pyruvate metabolism (xla00620)' (6 proteins, 
P=9.86xl0~ 3 ), and 'arginine and proline metabolism (xla00330)' 
(6 proteins, P=9.48xl0~ 3 ), were cold-exposure-specific. Minimal 
difference was observed between the control and cold-exposure 
conditions, possibly because of the small number ofX. laevis proteins 
annotated with GO and KEGG pathways. We then converted the 
identified proteins to their homologous human proteins for analysis 
(Fig. 4A) and found that 14 and 19 GO:BP terms were significantly 
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Fig. 5. Distribution of fold change in protein 
abundance. The bar chart shows protein ratios between 
the cold-exposure and control conditions for all 81 
relatively quantified proteins. 
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enriched in the control and cold-exposed X. laevis liver proteomes, 
respectively (Fig. 4B). Twelve of these terms were shared between 
control and cold exposure, and two and seven terms were control- and 
cold-exposure-specific, respectively (Fig. 4B; supplementary 
material Table S3). Five pathways were commonly associated with 
both control and cold-exposed X. laevis liver proteomes (Fig. 4C; 
supplementary material Table S4). 

Differential protein abundance under cold exposure 
To identify differentially expressed proteins associated with cold 
exposure, non-label quantification was performed (Fig. 5). The 
proteins found to be differentially expressed were categorized 
into four groups. Fourteen proteins were upregulated by cold 
exposure (fold change > 1 .25) (group 1; Table 1), and 13 proteins 
were detected only in the cold-exposure condition (i.e. 
considered to be newly induced) (group 2; Table 2). Nine 
proteins were downregulated under cold exposure (fold change 
<0.8) (group 3; Table 3), and the levels of 22 proteins were 
below the detection limit (i.e. detected only in the control 
condition) (group 4; Table 4). Generally, proteins in groups 1 and 

2 (total, 27) can all be viewed as upregulated and those in groups 

3 and 4 (total, 31) can be viewed as downregulated under cold 
exposure. We then analysed the GO enrichment of these 
upregulated and downregulated proteins by using the DAVID 



program with X. laevis proteins and candidate human 
homologues. Neither enriched GO:BP terms nor associated 
KEGG pathways were identified, presumably because of the 
small number of proteins. 

Validation of proteome data and GO and pathway analysis by 
metabolite measurements 

In the GO analysis, the 'tyrosine catabolic process' and '1- 
phenylalanine catabolic process' were significantly enriched only 
under cold exposure (Fig. 4B). In addition, levels of 
fumarylacetoacetase, which is necessary for metabolism of 
phenylalanine and tyrosine, tended to increase in cold exposure 
(Table 1). Generally, the liver metabolizes aromatic amino acids, 
including phenylalanine and tyrosine but not tryptophan (Dejong et 
al., 2007). Therefore, levels of these free amino acid levels in the 
liver were compared between the cold-exposure and control groups 
(Fig. 7A). Phenylalanine and tyrosine levels significantly decreased 
in cold exposure, whereas the tryptophan level did not change. 

Fig. 6 shows cold-exposure-induced changes in protein 
abundance associated with carbohydrate metabolic pathways. 
Glycogen phosphorylase (PYGL), which catalyses the 
breakdown of glycogen, was increased under cold exposure 
(Fig. 6). Levels of 1,4-alpha-glucan branching enzyme (GBE) 
and UTP-glucose-1 -phosphate uridylyltransferase (UDP-glucose 



Table 2. Newly induced proteins in cold-exposed X. laevis liver (group 2). 

Human homologous protein 







Protein 


Peptide 








Accession number 


Protein name 


score" 


content" 


RefSeq protein ID* 


Protein name 


Function 


gi 147905834 


lysosomal thioesterase 
PPT2-B precursor 


58 


11 


NP_005146 1 " 


lysosomal thioesterase 
PPT2 isoform 
a precursor 


protein-lipid modification 


gi 160420189 


dopa decarboxylase 
(aromatic L-amino acid 
decarboxylase) 


32 


6 


NP_00 1076440 


aromatic-L-amino-acid 
decarboxylase isoform 1 


amino acid metabolism 


gi 147898691 


DNA-dependent protein 
kinase catalytic subunit 


33 


6 


NP_008835 1 " 


DNA-dependent protein 
kinase catalytic subunit 


DNA repair, protein 
phosphorylation, cell cycle 










NP_005073 t 


isoform 1 


control 


gi 125858908 


unknown (protein for 
IMAGE:8550378) 


30 


5 


E3 ubiquitin/ISG15 
ligase TPJM25 


proteolysis 


gi 147902854 


zinc finger protein 507 


31 


5 


NP_055725 t 


zinc finger protein 507 


transcriptional regulation 


gi 147901600 


liver glycogen 
phosphorylase 


35 


4 


NP_002854 


glycogen phosphorylase, 
liver form isoform 1 


glycogen metabolism 


gi 68534041 


LOC733291 protein 


55 


4 


NP_006301 


puromycin-sensitive 
aminopeptidase 


proteolysis 


gi 147907415 


ropporin- 1 -like protein 


28 


3 


NP_001188395 1 ' 


ropporin- 1 -like protein 


unknown 


gi 148225833 


glycerol-3-phosphate 
dehydrogenase 1 (soluble) 


54 


3 


NP_005267 


glycerol-3 -phosphate 
dehydrogenase [NAD+], 
cytoplasmic isoform 1 


phospholipid metabolism 


gi 3745759 


histone H4 


75 


3 


NP_778224 t 


histone H4 


chromatin packaging and 
remodeling 


gi 147905626 


alcohol dehydrogenase 1C 
(class I), gamma 


35 


2 


NP_000660 


alcohol dehydrogenase 
1C 


carbohydrate metabolism 




polypeptide 






NP_002826 1 " 






gi 37791449 


protein tyrosine 
phosphatase PTP-PEST 


37 


2 


tyrosine-protein 
phosphatase 
non-receptor type 12 
isoform 1 


protein phosphorylation, 
signal transduction, cell 
adhesion, protein targeting, 
cell communication, cell 
motility 


gi 49115587 


LOC443650 protein 


32 


2 


NP_055807 


disheveled-associated 


cell motility 



activator of 
morphogenesis 1 



a Protein scores and peptide contents in MASCOT MS/MS ion search. Maximum values for triplicate data are shown. 
b RefSeq IDs of human homologues were obtained from NCBI HomoloGene database or by alignment with the NCBI Reference 
using the BLASTp programme ('). 



Sequence (RefSeq) database 
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Table 3. Downregulated proteins in cold-exposed X. laevis liver (group 3). 



Control 



Cold exposure 



Human homologous protein 



Accession 




Protein 


Peptide 


Protein 


Peptide 


Fold 


Refseq 






number 


Protein name 


score 


content 


score 


content 


change 


protein ID 


Protein name 


Function 


gi 33585701 


MGC69098 protein 


49 


3 


43 


2 


0.631 


NP_005882 


aceryl-CoA 
acetyltransferase 
cytosolic 


protein acetylation 


gi 148236249 


prosaposin precursor 


74 


6 


56 


3 


0.677 


NP_002769 1 ' 


proactivator 
polypeptide 

preproprotein 


lipid and fatty acid 
transport 


gi 147902842 


annexin A13 


75 


4 


75 


3 


0.693 


NP_001003954 


annexin A13 

1 jU1\_71 111 U 


cell differentiation 


gi 148233183 


elongation factor 1- 
alpha, somatic form 


159 


8 


127 


7 


0.728 


NP_001393 + 


elongation factor 
1 -alpha 1 


translational 
regulation 


gi 148226440 


MGC82879 protein 


547 


29 


445 


24 


0.730 


NP_001144 


annexin A4 


lipid, fatty acid and 
steroid metabolism 


gi 113571 


serum albumin B 


94 


6 


103 


5 


0.750 


NP_000468 1 ' 


serum albumin 
preproprotein 


transport 


gi 148237649 


sorbitol dehydrogenase 


76 


8 


41 


4 


0.750 


NP_003095 


sorbitol 
dehydrogenase 


carbohydrate 
metabolism 


gi 16332351 


glutathione 
S-transferase mu 2 


466 


21 


360 


16 


0.750 


NP_000839 


glutathione 
S-transferase Mu 2 
isoform 1 


detoxification 


gi 148222492 


uncharacterized protein 
LOC495316 


203 


11 


96 


8 


0.761 


NP_001087 


ATP-citrate 
synthase isoform 1 


lipid biosynthesis 



Columns are as described in Table l's footnote. 



pyrophosphorylase, UGPase), which participate in glycogen 
synthesis from glucose, were decreased (Fig. 6). Consequently, 
it is likely that synthesis of glycogen was downregulated and 
breakdown of glycogen was upregulated by cold exposure. Thus, 
we compared liver glycogen levels between the cold-exposure 
and control groups. The liver glycogen level tended to decrease 
to 70% of the control level upon cold exposure (Fig. 7B). 

Glycogen is broken down into glucose- 1 -phosphate (G1P) by 
PYGL and is subsequently converted to glucose 6-phosphate 
(G6P) by phosphoglucomutase (PGM). The resulting G6P is used 
in several metabolic pathways: (1) it is used in glucose 
production to supply energy for the body; (2) it fuels 
downstream glycolysis to produce energy in the form of ATP 
and NADH; (3) it is the starting substrate for the oxidative arm of 
the pentose phosphate pathway (PPP); (4) it may be used for 
biosynthesis of other metabolites. We examined whether liver 
glycogen was broken down for these pathways as follows: 

(1) Glucose production: increased PYGL and decreased 
hepatic glycogen were observed (Figs 6, 7B). PYGL is 
upregulated in winter freeze-tolerant wood frogs (Rana 
sylvatica), which produce glucose from glycogen as a 
cryoprotectant in freezing (Kiss et al., 2011). When hepatic 
glucose production is increased, glucose levels in both liver and 
plasma are elevated because the glucose moves to the 
bloodstream across the hepatocyte membrane, mainly through 
facilitated diffusion through the glucose transporter (Nordlie et 
al., 1999). Therefore, we compared liver and plasma glucose 
levels between the control and cold-exposure conditions. The 
liver glucose level was significantly decreased to half of the level 
in the control under the cold-exposure condition (Fig. 8A). In 
contrast, the plasma glucose level significantly increased in the 
cold-exposure condition (approximately 3.5 fold) (Fig. 8B). 

(2) Glycolysis: the level of pyruvate kinase (PK), which catalyzes 
the last step in the process of glycolysis by metabolising glucose to 



pyruvate with a net gain of two ATP and two NADH molecules per 
glucose molecule, was decreased (Fig. 6). 

(3) PPP: the 'pentose phosphate pathway' was commonly 
associated with the liver proteome in the control condition and 
under cold exposure (Fig. 4C). In addition, the 'pentose 
biosynthetic process' was significantly enriched only in 
under exposure (Fig. 4B). These results suggest that G6P 
derived from glycogen may be used for PPP. The oxidative 
phase of PPP generates NADPH from NADP (Wamelink et al., 
2008). \nX. laevis, cold exposure causes hepatic iron accumulation 
as a result of heme degradation following erythrocyte destruction 
(Maekawa et al., 2012). Iron accumulation is associated with free 
radical production and elevates oxidative stress (Bacon and 
Britton, 1990). Cellular NADPH is important for tolerance to 
oxidative stress and maintenance of cellular redox homeostasis. 
Therefore, hepatic NADP and NADPH levels were compared 
between the control and cold-exposure conditions and found to be 
unaffected by cold exposure (Fig. 8C,D). The hepatic NADPH/ 
NADP ratio was also not changed (data not shown). 

(4) Biosynthesis of other metabolites: in rainbow smelt 
(Qsmerus mordax), glycerol accumulation in the plasma is 
induced by low temperature (Driedzic et al., 2006). Some freeze- 
tolerant terrestrial anurans such as grey tree frogs {Hyla 
versicolor: and Hyla chrysoscelis) also produce copious 
quantities of glycerol as the cryoprotectant (Schmid, 1982; 
Irwin and Lee, 2003). These species use liver glycogen as the 
carbon source for glycerol synthesis. In this synthetic pathway, 
dihydroxy acetone phosphate (DHAP) is converted to glycerol 3- 
phosphate (G3P) and subsequently to glycerol via reactions 
catalysed by glycerol-3-phosphate dehydrogenase (GPDH) and 
glycerol-3-phosphatase (G3Pase), respectively. In our proteomic 
analysis, GPDH levels increased under cold exposure (Fig. 6), 
which suggests that glycogen may be used for glycerol 
production. To test this possibility, glycerol levels in the liver 
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Table 4. Proteins below detection limit in cold-exposed X. laevis liver (group 4). 



Human homologous protein 







Protein 


Peptide 


RefSeq protein 






Accession number 


Protein name 


score 


content 


ID 


Protein name 


Function 


m 14K737QQ7 


pnnlncp 1 n 

CTlUJltl&C 1 <x 


105 


9 


NP 001419 

111 \J\J liLy 


ri Inn n-pnnl n cp icnfrvrm 1 
ui una rjiiuidat 1 juiuini 1 


q1 vpnl'vci <; 


gi 148235435 


uncharacterized protein 


99 


6 


NP_006750 


UTP— glucose- 1 -phosphate 


polysaccharide metabolism 




LOC398814 








uridylyltransferase isoform a 




gi 148237900 


adenosylhomocysteinase B 


82 


3 


NP 000678 


adenosylhomocysteinase isoform 1 


purine metabolism 


gi 50415517 


heat shock 70 kDa 


78 


3 


NP_006588 


heat shock cognate 71 kDa protein 


protein folding, stress 




protein 








isoform 1 


response 


gi 50417653 


LOC397850 protein 


75 


3 


NP_005338 


78 kDa glucose-regulated protein 


protein folding, stress 












precursor 


response 


gi 288557262 


heat shock cognate 70.11 


71 


4 


NP 068814 1 


heat shock-related 70 kDa protein 2 


protein folding, stress 




rvrntpm 










rp cnflT15P 
ICTSIJVJIIjC 


gi 159155766 


LOC780759 protein 


69 


5 


N- 


tran sketol ase 


pentose-phosphate shunt, 










P_001 128527 




vitamin metabolism 


gi 147899332 


glucan ( 1 ,4-alpha-), 


63 


6 


NP_000149 


1 ,4-alpha-glucan-branching enzyme 


glycogen metabolism 




branching enzyme 1 












gi 11493740 


DNA polymerase 


50 


5 


NP 006222 1 


DNA. polymerase epsilon catalytic 


DNA repair, cell cycle 




epsilon 






NP 000272 f 


subunit A 




gi 18202614 


ptenn-4-alpha-carbinolaiTiine 


48 


2 


pterin-4-alpha-carbinolamine 


mRNA transcription, pterin 




npnvni'atasp 








L1C11 V Ul£l Lao^ 


tnptnnn 1 cm 


gi 64659 


plnncmtinn friptnr 1 -nlnlm 

ClUlliiClllVJll 1&UIVJ1 1 tlllJlla 


47 


3 


NP 001393 


plnriCTnti nn fnctrvr 1 -nlnrifi 1 
CHJllliaLlUll lutLUl 1 allJlla. l 


trnncl n tinnnl vpctii 1 ntinn 

11 ull^ld LIUllul 1L_1.1K111VM1 


gi 148233056 


pyruvate kinase liver and 


46 


2 


NP 002645 


pyruvate kinase isozymes IVI1/M2 


glycolysis 




RBC 








i cnffi rm n 
lOVJlVJllll a 




gi 147902366 


cathepsin C precursor 


39 


4 


N- 


dipeptidyl peptidase 1 isoform c 


proteolysis 










P_001107645 t 


precursor 




gi 147905009 


pericentriolar material 1 


38 


9 


NPjioeiss 1 " 


pericentriolar material 1 


centrosome assembly 




protein 






NP_002412 t 


protein 




gi 71051388 


matrix metalloproteinase-18 


38 


3 


interstitial collagenase 


proteolysis 




precursor 








isoform 1 preproprotein 




pi 147907076 

IT / 7 V7Z. u/. KJ 


npfnYit'pnnviti f\ 

UC1VJA.11 CLIVJA.111 \J 


38 


3 


NP 004896 


nprriYivprlfiYiTi-fi 

jJCl UAlltUUAlll u 


iinmiinih/ nnn npfpnsp 

1111111LL111L V clllU UClCll^VT, 














antioxidation and free radical 














removal 


pi 147906799 

iil It / y \J\J j y y 


Iptnnl cunnt lan;ap nnmnlAO 1 
l&lllal glullL lal Varj 11U111U 1U ti l 


36 


5 


NP 004 1 3 1 

INI J 1 


lptnnli'?^ CTinnt lurvnp nrnt"Pin 

1 L- 1 1 Kl l\ 1 1 u 1 1 L lal Va^ 111 U 1 L 1 1 1 


TPCTll 1 fltPfl PYflPlftfl*;! S 
IL—LIKlld.! CAUL VIUdIB, 












homolog 1 


neurogenesis, ectoderm 














development 


gi 148222987 


glucoside xylosyltransferase 


35 


1 


N- 


glucoside xylosyltransferase 


O-glycan processing 




2 precursor 






P_001073862 t 


2 precursor 




gi 147906817 


tumor necrosis factor, 


33 


4 


NP_006281 t 


tumor necrosis factor 


anti-apoptosis, inflammatory 




alpha-induced protein 3 








alpha-induced protein 3 


response, ubiquitin 














conjugation pathway 


gi 148229521 


protein CIP2A homolog 


32 


3 


NP_065941 


protein CIP2A 


unknown 


gi 147904142 


topoisomerase (DNA) II 


31 


15 


NP_008958 


DNA topoisomerase 


DNA replication, cell cycle, 




binding protein 1 








2-binding protein 1 


chromatin packaging and 














remodeling 


gi 51703567 


MGC84331 protein 


31 


7 


NP_078857 


protein FAM184A isoform 1 


unknown 



Columns are as described in Table 2's footnote. 



and plasma were compared between the control and cold- 
exposure conditions. Glycerol levels in both the liver and plasma 
were not changed by cold exposure (Fig. 8E,F). 

Discussion 

Here, using a proteomics approach, we demonstrated the protein 
abundance profile of X. laevis liver after cold exposure. One 
hundred and forty-five proteins were identified from the soluble 
fraction of whole liver extract by trypsin digestion followed by 
LC-MS/MS without prior fractionation of proteins; 58 of these 
proteins were differentially expressed proteins (upregulated, 27; 
downregulated, 31). No GO:BP terms, however, were significantly 
enriched in these protein lists, presumably due to the small number 
of proteins. In recent years, global proteomics approaches have 
been widely used to characterize a number of tissue proteomes, 
including liver (Falcon-Perez et al., 2010). Pre-fractionation of 
protein samples by one-dimensional electrophoresis and/or liquid 



chromatography has enabled the identification of more than 7000 
proteins in the mouse liver (Shi et al., 2007; Lai et al., 2008). By 
the application of these methods, the number of identified proteins 
will be increased and enable GO and pathway analysis. 

We have shown that the number of mature erythrocytes and 
accumulation of iron catalyzed from heme decrease in the liver 
within 24 hours after cold exposure (5°C), which reflects the 
enhanced destruction of erythrocytes in the liver (Maekawa et al., 
2012). The enhanced hepatic destruction of erythrocytes is 
thought to increase the protein levels of haemoglobin subunits in 
the liver. Consistent with a previous finding, the protein levels of 
haemoglobin subunits in the liver tended to increase within 
24 hours after cold exposure in our proteomic analysis (Table 1), 
which suggests that the results of the LC-MS/MS analysis in our 
study reflect physiological responses as changes in protein 
abundance. The peripheral erythrocyte count remains low during 
cold exposure (5 days) because of the hepatic confinement of 
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Fig. 6. Cold-exposure-indueed changes in protein abundance associated with carbohydrate metabolism. Substrates and enzymes are as follows: 
(glycolysis/glyconeogenesis) G6P, glucose 6-phosphate; F6P, fructose 6-phosphate; FBP, fructose 1 ,6-bisphosphate; GAP, glyceraldehyde 3-phosphate; 
DHAP, dihydroxyacetone phosphate; 1,3BPG, 1,3-bisphosphoglyceric acid; 3PGA, 3-phosphoglyceric acid; 2PGA, 2-phosphoglyceric acid; PEP, 
Phosphoenolpyruvate; GCK, glucokinase; G6Pase, glucose-6-phaosphatase; GPI, glucose-6-phosphate isomerase; PFK, phosphofructokinase; FBPase, 
fructose- 1 ,6-bisphosphatase; FBA, fructose 1 ,6-bisphosphate aldorase; TPI, triosephosphate isomerase; GAPDH, glyceraldehyde-3-phosphate dehydrogenase; 
PGK, phosphoglyceric acid kinase; ENO, enolase; PEPCK, phosphoenolpyruvate carboxykinase; PK, pyruvate kinase; MDH, malate dehydrogenase; LDH, 
lactate dehydrogenase; ACL, ATP-citrate synthase (ATP-citrate lyase); ACS, acetyl-CoA synthase; ALDH, aldehyde dehydrogenase; ADH, alcohol 
dehydrogenase; (pentose phosphate pathway) 6PG, 6-phosphogluconate; Ro5P, ribose 5-phosphate; Ru5P, ribulose 5-phosphate; Xu5P, xylulose 5-phosphate; 
S7P, sedoheptulose 7-phosphate; E4P, erythrose 4-phosphate; G6PDH, glucose-6-phosphate dehydrogenase; 6PGDH, 6-phosphogluconate dehydrogenase; 
TKET, transketolase; TALD, transaldolase; (glycogen metabolism) G1P, glucose 1-phosphate; UDPG, uridine diphosphate glucose; PGM, phosphoglucomutase; 
UGPase, UDP-glucose pyrophosphorylase; GYS, glycogen synthase; PYGL, glycogen phosphorylase; GBE, glycogen branching enzyme; GDE, 
glycogen debranching enzyme; (others) Glycerol-3-P, glycerol 3- phosphate; F1P, fructose 1-phosphate; GK, glycerol kinase; G3Pase, glycerol-3-phosphatase; 
GPDH, glycerol-3-phosphate dehydrogenase; FK, fructokinase; ALDB, fructose-bisphosphate aldolase B; GAK, glyceraldehyde kinase; GRHPR, 
glyoxylate reductase/hydroxypyruvate reductase; TMABADH, 4-trimethylaminobutyraldehyde dehydrogenase; ACAT, acetyl-CoA acetyltransferase; 
mMDH, mitochondrial MDH; IDH, isocitrate dehydrogenase. Modified from portions of KEGG pathway map for 'glycolysis/gluconeogenesis' (00010) 
and 'Pentose phosphate pathway' (00030). 
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Fig. 7. Comparisons of free aromatic amino acid levels and glycogen levels 
in the liver between the control and cold-exposure conditions. (A) Free 
amino acids in the liver («=3). (B) Glycogen in the liver («=4). Each bar 
represents the s.e.m. *P<0.05 by Student's Mest. 

newly produced erythrocytes (Maekawa et al., 2012). Some of 
the differentially expressed liver proteins may play a role in this 
phenomenon. Long-term exposure to 10°C (153 days) also causes 
chronic pancytopenia in X. laevis (Maekawa et al., 2012). Our 
proteomics approach is probably applicable to the investigation 
of the response to cold over a longer period. 

Consistent with the increase in fumarylacetoacetase and cold- 
exposure-specific GO:BP terms ('tyrosine catabolic process' and 
'1-phenylalanine catabolic process'), free phenylalanine and 
tyrosine levels in the liver decreased (Table 1, Figs 4B, 5 A). 
This finding suggests that amino acids were catabolized for 
energy, ketogenesis, and/or gluconeogenesis, and it shows the 
utility of a proteomics approach including GO analysis inX. laevis. 

Increased PYGL and decreased hepatic glycogen were 
observed (Table 2, Figs 6, 7B). PYGL is increased in winter 
freeze-tolerant wood frogs (Rana sylvatica) that produce glucose 
as a cryoprotectant from hepatic glycogen upon freezing (Kiss et 
al., 2011). Rainbow smelt (Osmerus mordax) and some freeze- 
tolerant terrestrial anurans such as grey tree frogs (Hyla 
versicolor and Hyla chrysoscelis) produce glycerol as a 
cryoprotectant from hepatic glycogen in winter (Driedzic et al., 
2006; Schmid, 1982; Irwin and Lee, 2003). Therefore, we 
considered that glycogen breakdown in the liver of X. laevis 
results in production of glucose or glycerol. The glucose level in 
the liver, however, was decreased and liver and plasma glycerol 
levels were not changed (Fig. 7A,E,F). Therefore, glycogen was 
utilized in neither glucose nor glycerol production. Levels of PK, 
which participates in glycolysis, were decreased upon cold 
exposure, which suggests that glycolysis was inhibited by cold 
exposure and that G6P derived from glycogen was not used for 
energy production. Therefore, it is possible that the breakdown of 




Control Cold 
exposure 



Control 



Cold 
exposure 



Fig. 8. Comparisons of glucose, glycerol, NADP, and NADPH levels in the 
liver and/or plasma between the control and cold-exposure conditions. 

(A,B) Glucose levels in the liver («=3) and plasma (« = 5). (C,D) NADP and 
NADPH levels in the liver (n=4). (E,F) Glycerol levels in the liver («=4) and 
plasma («=4). Each bar represents the s.e.m. *P<0.05 by Student's Mest. 

glycogen is utilized in other pathways. Even though cold-exposed 
X. laevis presents hepatic iron accumulation that causes oxidative 
stress, levels of NADPH, which is consumed for protection 
against oxidative damage, were not decreased (Fig. 7C). This 
finding suggests that G6P derived from glycogen may be utilized 
for NADPH production by PPP for protection against oxidative 
damage (Fig. 9). Considering that cold-induced apoptosis in 
cultured rat hepatocytes and liver endothelial cells is mediated by 
reactive oxygen species (Rauen et al., 1999), provision of 
NADPH for protection from oxidative stress in the liver of cold- 
exposed X. laevis may provide a mechanism for coping with cold 
in freeze-intolerant ectothennic vertebrates. 

The plasma glucose level was increased although the hepatic 
glucose level decreased (Fig. 7B). Considering that the liver is 
the only organ that actually releases glucose into the systemic 
circulation, our findings may reflect a decrease in circulating 
erythrocytes that consume glucose and decreased tissue glucose 
consumption resulting from decreased metabolism under low 
body temperature. 

Levels of the endoplasmic reticulum (ER) molecular 
chaperone, 78-kDa glucose-regulated protein (GRP78, also 
known as BiP or Hspa5), were decreased under hypothermia in 
the liver of X. laevis (Table 4). GRP78 was found to be similarly 
decreased in freeze-tolerant wood frogs in the winter (Rana 
sylvatica) and in hypothermic rat livers (Oda et al., 2012; Kiss et 
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Fig. 9. Schematic models of the early response to cold exposure in the liver 

of X. laevis. In cold exposure, hepatic glycogen and glucose are thought to be 
used in the pentose phosphate pathway for NADPH supply rather than in energy 
production through glycolysis. This mechanism suppresses oxidative stress 
derived from iron accumulation caused by hepatic erythrocyte degradation. 
Abbreviations for the substrate and enzymes are shown in Fig. 6. 

al., 2011). In contrast, the content of GRP78 was found to be 
increased among the liver proteins of hibernating mammals 
(Epperson et al., 2004; Epperson et al., 2010). GRP is induced 
during ER stress derived from the disruption of calcium 
homeostasis and accumulation of unfolded proteins in the ER, 
and it protects the host cell against ER stress-induced cell death 
(Rao et al., 2002). Increased GRP expression is unique to 
hibernating animals and appears to be important. 

This study presents the first system-wide screen of hepatic 
proteins in X. laevis. Low environmental temperature (i.e. low 
body temperature in X. laevis) induced modulation of the hepatic 
proteome, particularly in pathways associated with carbohydrate 
metabolism and with other functions. Regarding carbohydrate 
metabolism, this study suggests that utilization of hepatic 
glycogen for NADPH supply is associated with tolerance to 
cold-induced oxidative stress. 

This study not only provides insight into the mechanism of 
tolerance to low body temperature but also demonstrates that 
proteomic analysis of X. laevis is applicable despite the lack of 
whole-genomic information because of the accumulation of 
expressed sequence tags (ESTs) in public databases. In X. laevis, 
the liver is a major site for both metabolism and haematopoiesis. 
Therefore, the modulation of the hepatic proteome related to 
metabolic pathways under low body temperature occurs at a 
haematopoietic site. This study provides a model for 
investigation of not only coping mechanisms for cold but also 
interactions between haematopoietic and metabolic systems. 
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